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ABSTRACT 

The effects of crossflow on the growth rate of inviscid Gortler vortices in a hypersonic 
boundary layer with pressure gradient are studied in this paper. Attention is focussed on the 
inviscid mode trapped in the temperature adjustment layer; this mode has greater growth 
rate than any other mode. The eigenvalue problem which governs the relationship between 
the growth rate, the crossflow amplitude and the wavenumber is solved numerically, and the 
results are then used to clarify the effects of crossflow on the growth rate of inviscid Gortler 
vortices. It is shown that crossflow effects on Gortler vortices are fundamentally different for 
incompressible and hypersonic flows. The neutral mode eigenvalue problem is found to have 
an exact solution, and as a by-product, we have also found the exact solution to a neutral 
mode eigenvalue problem which was formulated, but unsolved before, by Bassom and Hall 
(1991). 
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1 Introduction 


This paper continues our studies on the Gortler instability mechanism in hypersonic bound- 
ary layers. For a discussion of the motivation for such studies and background information, 
the reader is referred to our previous papers, Hall and Fu (1989) and Fu, Hall and Blackaby 
(1992). Both are concerned with linear stability properties of Gortler vortices, but the first 
studies a Chapman’s law fluid and the second is devoted to a Sutherland’s law fluid. A com- 
mon feature between the results for these two types of fluids is that the most rapidly growing 
Gortler vortices in a hypersonic boundary layer are those trapped in the temperature adjust- 
ment layer sitting at the edge of the boundary layer over which the basic state temperature 
decays from its O ( M 2 j values near the wall to its free stream value, here M is the free stream 
Mach number. The major differences when these two different viscosity laws are used are as 
follows. Firstly, using Chapman’s viscosity law predicts an 0(M 2 ) Gortler number, whilst 
using Sutherland’s viscosity law predicts a smaller 0(M 3 ^ 2 ) Gortler number. Secondly, the 
evolution of Gortler vortices in Chapman fluids is less strongly affected by nonparallel effects 
than those in Sutherland fluids. This is due to the fact that when Chapman’s law is used, the 
temperature adjustment layer is logarithmically thin whilst when Sutherland’s law is used, 
the adjustment layer is much thicker. Certainly, for a hypersonic boundary layer across 
which the temperature varies significantly, it is more appropriate to use Sutherland’s law 
to model the viscosity of the fluids. Thus, in our later papers, Fu and Hall (1992a, 1992b) 
which study nonlinear development and secondary instabilities, Sutherland’s viscosity law 
has been consistently used. In Fu and Hall (1992a), we have shown how Gortler vortices may 
grow in the neighbourhood of the neutrally stable position to form a large amplitude vortex 
structure which consists of a core region bounded by two viscous transition layers. Once 
such a structure has been established, various secondary instabilities in the form of traveling 
waves may be triggered off. There are at least two such secondary instabilities which have 
been observed in experiments. The first is called the wavy vortex mode of instability which 
makes the vortex boundaries wavy and which is located in the two viscous transition layers. 
This secondary instability has been studied in Fu and Hall (1992a). The other type is an 
inviscid instability which can exist because of the inflexional nature of the steady state. In 
Fu and Hall (1992b), we have shown how the change in the basic flow caused by the large 
mean flow correction of Gortler vortices would affect the growth rate of Rayleigh instability. 

All of our previous studies in the hypersonic context have been concerned with problems 
which arise when Gortler vortices occur i n two-dim ensional boundary layers, but in many 
practical situations, the basic boundary layer may be three-dimensional. Boundary layer 
flows over a three-dimensional obstacle or over a turbine blade are such examples. A more 
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significant example is the boundary layer flow over the two areas of concave curvature on the 
lower side of a modern laminar flow airfoil which is fully three-dimensional when the wing 
is swept. Obviously, in order to effectively control or delay transition, the role played by a 
crossflow component in the basic flow should be fully understood. 

This paper may be considered as extending the Hall (1985) and Bassom and Hall (1991) 
studies on crossflow effects in incompressible flows to hypersonic flows. It is also related to 
the recent work by Dando (1992) on crossflow effects in compressible flows. All of these 
studies have shown that crossflow has a stabilizing effect on Gortler vortices, and more 
importantly, Gortler vortices may be completely stabilized by the three-dimensional nature 
of the basic state. One objective of the present paper is to give a quantitative description of 
the stabilizing effect of a crossflow on Gortler vortices in hypersonic boundary layers. 

As has been observed before by Hall (1985) for incompressible flows, the case when the 
crossflow velocity component is a constant multiplier of the streamwise velocity component 
is degenerate in the sense that by a suitable transformation, the problem can be reduced to 
a Gortler problem corresponding to a two-dimensional basic state. This degeneracy happens 
when there is no pressure gradient in the inviscid flow. Thus in order to reveal the genuine 
implications of having a crossflow, we avoid this degeneracy by assuming in this paper 
that there exists a pressure gradient in the inviscid flow. However, as has been shown 
by Stewartson (1964), when there is a pressure gradient in the boundary layer, analytical 
solutions to the boundary layer equations are possible only if we make three assumptions, 
one of which is that the viscosity should be given by Chapman’s law. In order to make the 
necessary analytical progress, we shall assume that this is the case. Thus the present paper 
is more relevant to Hall and Fu (1989). 

An important property of a hypersonic boundary layer is that the chordwise acceleration 
of a fixed fluid particle in the basic state is asymptotically large and in the opposite direction 
to the centrifugal acceleration due to the curvature of the wall. This negative chordwise 
acceleration has been interpreted in our previous papers as effectively producing a negative 
curvature (or equivalently a negative Gortler number). Thus the curvature of the wall should 
be large enough to overcome this negative curvature due to the basic state if Gortler vortices 
are to exist at all. In this paper, we assume that the Gortler number is of the same order 
as that of the effective negative Gortler number due to the basic state curvature. Then the 
large resultant Gortler number drives a large growth rate and as a result, the eigenvalue 
problem determining the growth rate is essentially inviscid. The main aim of this paper is 
to determine how the existence of a crossflow would modify the growth rate of these inviscid 
Gortler vortices. 

The rest of this paper is organized as follows. As a preparation to later analysis, we first 
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discuss in the next section the solutions to the three-dimensional basic state when there is 
a pressure gradient. In section 3 we formulate the general linear Gortler instability problem 
which is then specialized in section 4 to the inviscid problem mentioned above. Sections 
4 and 5 are devoted to the non-neutral and neutral solutions respectively, and in the final 
section we discuss our results and make some conclusions. 


2 The basic state 


Consider a hypersonic boundary layer over a rigid wall of variable curvature { \ ( A)k{x* j L), 
where L is a typical streamwise length scale and A is a lengthscale characterizing the radius 
of curvature of the wall. We choose a curvilinear coordinate system (x*,y*,z*) with x* 
measuring distance along the wall, y* perpendicular to the wall and z* in the spanwise 
direction. To fix ideas, we may assume that the rigid wall is the surface of an infinite 
cylinder; the z-axis is then aligned with a generator of the cylinder (see Fig.l). The velocity 
components are denoted by (u*, v * , w*) and density, temperature and viscosity by p*, T* and 
fi* respectively. The free stream values of these quantities will be signified by a subscript oo. 
The Reynolds number R, the curvature parameter A and the Gortler number G are defined 

by 

U cc L Poc A L Q _ 2#l/2 A 


R = 


A = 






( 2 . 1 ) 


We assume that the Reynolds number is large, while A is taken to be small. In the context 
of incompressible flows, the characteristic values of G for instability are 0(1); whilst in 
hypersonic flows G is 0(M 2 ) when Chapman’s viscosity law is used (see Hall and Fu (1989)), 
and is 0(M 3 / 2 ) when Sutherland’s viscosity law is used (seeFu, Hall and Blackaby (1992)). In 
the following analysis, coordinates (x* ,y* , z*) are scaled on (L, R~ X ^ 2 L, R~*/ 2 L), the velocity 
{u*,v*,w*) is scaled on (ix^, Rr^^u^, R~ l ^ 2 ul Q ) and other quantities such as p*,T m , and /x* 
are scaled on their free stream values with the only exception that the pressure p* is scaled 
on and the bulk viscosity A* is scaled on /x^. All dimensionless quantities will be 

denoted by the same letters without a superscript *. For an ideal gas without dissociation, 
the continuity, Navier-Stokes, energy equations and the equation of state are, to leading 
order, given by 

( 2 . 2 ) 
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Here terms of relative order R~ x have been neglected and we have used a mixed notation in 
which (vi,W 2 ) v 3) is identified with (u,v,w) and (xi,x 2 ,X3) with ( x,y,z ). Repeated suffices 
(3 signify summation from 1 to 3. The constants 7, M and <x are in turn the ratio of specific 
heats, the free stream Mach number and the Prandtl number defined by 
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where is the gas constant, k is the coefficient of heat conduction, and ^OO is 

the speed of sound in the free stream. In equations (2.3)— (2.6), the operator Dj Dt is the 
material derivative and is given by 

D d d d 

— = u— + v— + w—. 

Dt ox dy dz 

The basic state is given by 

(u, V, w) = (u(x,y), v(x,y), w(x,y)), T = T(x,y), 

p = p(x,y), p = ji(x,y), (2.8) 

where we have allowed for a non-zero crossflow component w(x,y). Such a three-dimensional 
boundary layer may easily be obtained if, for example, the angle of yaw of the cylinder A 
shown in Fig. 1 is non-zero. 

By substituting (2.8) into the governing equations (2.2) — (2.7), we find that the basic 
state satisfies 

A A 

(2.9) 
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'lM 2 P\ = /57\ (2.13) 

where P] is the inviscid solution for the pressure evaluated at y = 0. The equation obtained 
from the t/-momentum equation (2.4) is simply dp/dy = 0 which implies that p = p{x) = Pj . 
From the inviscid theory, we have the relation 


dPi TT dU\ 

= ~P\U\ 


(2.14) 


dx r ‘" dx’ 

where p\ and U\ have the same meaning as P\. If we introduce Ti, the inviscid solution for 
the temperature evaluated at y = 0, then equation (2.13) can also be written as 


pT = P\T\. 


(2.15) 


In this paper, we assume that dP\jdx ^ 0, i.e. there exists a pressure gradient in the 
inviscid flow. If dP x jdx — 0, then (2.10) and (2.11) together with the boundary conditions 
imply that w is a constant multiplier of u, in which case the original three-dimensional basic 
state can be made two-dimensional by rotating the coordinate system about the y-axis by 
an angle given by tan -1 (w/u), see Hall (1985). 

We now turn to the solution of the basic state equations. The reader is referred to the 
book by Stewartson (1964) for a detailed discussion of these equations. Here we shall just 
give an outline of the necessary procedures and the main results. 

The immediate consequence of including a pressure gradient is that similarity solutions 
are now only possible if we use Chapman’s viscosity law 

p = C(x)T, (2.16) 


assume that a = 1 and make the homenergic assumption. (We recall that analytical similar- 
ity solutions are possible without making any of the above assumptions if there is no pressure 
gradient). In this paper, we shall adopt these assumptions to make the necessary analytical 
progress. Then the first step in the solution is to determine the temperature T explicitly. To 
this end, we first eliminate dP\jdx from (2.10) and (2.12) to obtain the following equation: 
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where 


'=f+i( 7 -l)M 2 u 2 . 


(2.18) 


On comparing (2.17) with (2.11), we deduce that the general solution to (2.17) can be written 
as 

/ = —Aw + const., 

where A is an arbitrary constant to be fixed by boundary conditions. Hence 


T H — (7 — \)M 2 u 2 + Aw = const. 
2 

= r, + i(7 - 1)M 2 U 2 + AW U 


(2.19) 


where W\ is the inviscid crossflow velocity component evaluated at y = 0 and is given by 

Uryn sin A 


W l 


— R 1 / 2 tan A, 


since = Coo cos A (see Fig.l) and W\ has been scaled on /2 _1 ^ 2 u^ 0 . Guided by previous 
results for incompressible flows, we assume that the yaw angle A is of order R~ 1 ^ 2 (this is 
the order at which crossflow begins to affect the growth of Gortler vortices), so that is 
an order one constant under the large Reynolds number limit. We see that the the size of 
can be adjusted by varying A. 

By the homenergic assumption, the total energy (scaled by CpT^) in the inviscid flow is 
a constant and so 

r, + 5(7 - 1 )M 2 C , 2 = 1 + i (7 - 1)A/ 2 , (2.20) 

where the left and right hand sides are the scaled energy evaluated at y = 0 and y = oo, 
respectively. With the use of (2.20), (2.19) can be written as 


? + i ( 7 - lJA/V + A(w - Wi) = 1 + t( T - 1 )M 2 . 


( 2 . 21 ) 


If we denote the temperature at the wall by T w , then (2.21) together with u(0) = ie(0) = 0 
implies that 

T w - 1-|(7-1)M 2 


A 


Wi 


and so (2.21) becomes 
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r = i + §(7 -i)A/ 2 (i 


“ 2 ) + ( 1 -^){ r »- 1 -^- 1)M2 } 


( 2 . 22 ) 


(2.23) 


In the special case when the wall is insulated ( dT j dy = 0 at y = 0), it is easy to show with 
the aid of (2.23) that the wall temperature must necessarily be given by 


^Iwall insulated — ^ + nil 0^ T r , 


(2.24) 
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where T r is usually referred to as the recovery temperature. Once we know the expression 
for T, we can now solve (2. 9) -(2. 13) by introducing a variable Y by 


Y 


f y - 

= / pdy 

Jo 


(2.25) 


and expressing u and v in terms of a function tp(x,Y) as 

dtp _ l ( dtp 

U = W • * = ~pai + uY - ] - 

Then the continuity equation (2.9) is satisfied and (2.10), (2.11) become 

d 3 tp 

dY 3 ' 


(2.26) 


dPi 


d tp d 2 tp dtp d 2 tp ' 

WdYdY ~ lYxdY 2 = _ 7 M 2 P, dx 


+ C{x) 1 M 2 P \ 




(2.27) 


(2.28) 


Before we could look for similarity solutions for the above equations, we have to eliminate 
the appearance of u through T on the right hand side of (2.27). This can be achieved by 
introducing two new variable (X, Y) through 


X = J%M 2 P{x)^Tp7)C{x)dx, Y = yfp{x)Y. 

In terms of X and Y , equations (2.27) and (2.28) can be written as 

dtp d 2 tp dtp d 2 tp d 3 tp dVy 


dY dXdY dX dY 2 dY 3 

dtp dw dtp dw d 2 w 
dYdX ~ dXdY = W 2 ' 
after (2.15), (2.13), (2.14) and (2.20) have been used. Here 


+ v '7 p{' Jr[x ~ 1) }’ 


V t (X) = 


Ux 

VTt 


(2.29) 

(2.30) 

(2.31) 

(2.32) 


Equations (2.30) and (2.31) are to be solved subject to the following boundary conditions: 


* = ^ = 0, w = 0 at y = 0, 

dy 

?£->Vu — ► Wi as Y -» oo. 

dy 

The above boundary value problem admits similarity solutions of the form 

Y 


tp - a(X)f(r /), w = W\g(r}), rj = 


b(xy 


(2.33) 


(2.34) 
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if Vi is of the form 

Vi = UX m . 

By substituting (2.34) and (2.35) into (2.30) — (2.33), it is easy to show that 


>(*) = 


I 2U 

1 + m 


■X l +* 


b(X) = 


1/(1 -(- m) 




and the boundary value problem then reduces to the following simpler form: 

/"' + //" + fi{\ ~ f' 2 + (n - 1)(1 - 5 )} = 0, 

g" + fg' = o, 

/(0) = /'(0) = flr(0) = 0, ■ 

/'( oo.) = g (° °) = !> 

where the constant 0 and the cooling coefficient n in (2.37) are given by 

2m T,„ 


0 = 


1 + m '' 


n 


(2.35) 

(2.36) 

(2.37) 

(2.38) 

(2.39) 

(2.40) 

(2.41) 


On substituting (2.34a) into (2.26a), we obtain u = j) and insertion of this expres- 

sion and (2.34b) into (2.23) then yields the following expression for T: 


T = T,r + [1 + 1(7 - 1)M 2 ] {l - r + (» - 1)(1 - ?)} . 


(2.42) 


Our previous experience with the zero-pressure gradient problem shows that the stability 
properties of a hypersonic boundary layer depend very much on the properties of the basic 
state temperature distribution. When there is no pressure gradient, it is known from Hall 
and Fu (1989) that the boundary layer splits into two sublayers in the large Mach number 
limit: a wall layer where T = 0(M 2 ) and a logarithmically thin temperature adjustment 
layer where T = 0(1), and it is the latter sublayer which supports the most rapidly growing 
mode. In order to determine the precise boundary layer structure for the present non-zero 
pressure gradient problem, we first need to know the large r/ behaviour of /' and g. It can 
be shown from (2.37) and (2.38) that as rj — > oo, 

(n - l)c 2 


I'(V) = 1 + 


c i 


2(rj — cj) (g — v) 1+2 P 

5=1“ 
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C - 2 . e -§( r >~") 2 

Tj — U) 
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(2.43) 


where u>, c\ and c 2 are constants. On substituting (2.43) into (2.42), we find that as tj — » oo, 


T = Ti + NM 2 e-^~ w)2 • 


1 

(g-uY+W 


+ 


(2.44) 
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where 


N = ct{ 7 - 1 ). 


In order to determine the location and thickness of the temperature adjustment layer at 
which T = 0(1), it is convenient to let 


M 2 e -±( 7 ,- 0,) 2 . 


Solving this equation for (77 — to) gives 


T] — iO = Ml — 


= T± e y 


(77 -ujy+W N 


y + (1 + 2 /?) log Mi 
M, 


(2.45) 


where 

Mi = ^2 log M 2 . 

Equation ( 2 . 45 ) shows that the temperature adjustment layer is at 77 = 0(Mi) and is of 
0(1/Mi) thickness. In the adjustment layer, T is given by 


f = Ti(l+e s ), (2.46) 

so that T — > T\ as y — > —00 and matches with the wall layer solution (2.44) as y — > + 00 . 

In summary, the basic state under consideration corresponds to Vj = U\/\/T\ = UX m 
and has solution given by 

S = t7,m S = -^|^7? / V(l)o'(X)/(,) + C,r(,)(V; + ^5X)|, (2.47) 

IS = vr, <,(,), T = T,/' J + [1 + t(l - 1)M 2 ] {l — z' 2 + ^j,(l — s)} , (2.48) 

where the new variables X,Y,Y and rj are related to (x,y) by 

X = f 1 M 2 P x {x)Jr^)C(x)dx, Y = sJt&)Y, Y = j*pdy, 77 = (2.49) 


and it can be shown that the function Y x appearing in (2.47b) is given by 

n 7%, . i>» i 


= T(x,Y)d? -jj o 


dx 


(2.50) 


In section 4, we shall need to know the asymptotic expressions for these basic state quantities 
and their various derivatives in the temperature adjustment layer in order to determine the 
relative orders of the perturbation quantities. These expressions are listed in the Appendix 
at the end of this paper. 
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3 The linear Gortler instability problem 


In order to determine the evolutionary behaviour of Gortler vortices in the hypersonic bound- 
ary layer described above, we now superimpose a vortex structure on the basic state. Thus 
the total flow is written as 

u = u + Ue iaz , v = v + Ve iaz , w = w + We iaz , 


p=P 1+ —Pe iaz , T = T + 6e iaz , p = p + pe iaz . (3. 1 ) 

n 

Here U, V, etc. are functions of X and Y only. On substituting (3.1) into the governing 
equations (2.9)— (2.13) and linearizing, we obtain the following perturbation equations: 

f5 L (T) - - fin* + «») - f( L {0) + ia w0) 

*£*£—^2- 1)-' m 


p {L(U) + iawU + Uu x + Vuy) — ^ L(u)0 
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, -r- ,2-SU /, d 2 U _ 9U . fidW 

2 4 d ( 3 V 

+i.a • -p y W + p {L(V) + iawV} -f pv y V + a 2 fiV - Hz— 

OP dQ . „ .2 , 4 _ . 69 

+ Y ~ ^’Yx ~ 

~ if [t(») + jCkk 2 ] - ^ + (fru„) x \ 9 = 0, 


■ - rr . pdU __ TT . . pdV 4_ 2 

iap x U + ta ■ — — + iap y V -(- za • — — -pa W 

3 ox 6 ay 6 

6 ( 9W \ T/ . „ . 2., , 60 

+ dy \ ~dy' ) ~ pWyV ~ mF ~ ’ ^ Ux + Vy '° + pWy ~6y 


+ [f L i w ) + (P™y)yW ~ p[ L ( W ) + iawW] = 0, 

\i>T x - (7 - 1)V 2 ^-]C - 2(7 - 1 )M 2 - + pT y V 

+i .«i +•>.§»- 


10 



(3.6) 


- [tL(T) + (7 - i)M 2 fi(u y y + !(/;?„)„} e = o, 



Here fi = T ), y — dy/dT , and the operator L( ) is defined by 


(3.7) 


T , . _ 5 _ 9 

i() = “a; + V 


In these equations, variables x and y have not yet been transformed into X and Y using 
(2.49). This will be done at a later stage. 

So far we have not specified the size of the wavenumber a relative to the Mach number. 
We know from the analysis in §2 that the hypersonic boundary layer splits into two sublayers 
in the large Mach number limits: a wall layer of 0(M\) thickness where T — 0(M 2 ) and 
a temperature adjustment layer of 0(1) thickness sitting at the edge of the boundary layer 
where T — 0(1). We note from (2.49) that 


dy = 


KX) 

PiT ? /2 


Tdy. 


Hence in terms of the physical variable y , the wall layer and the temperature adjustment layer 
are of thickness of 0(M 2 M\) and 0(1/Mi), respectively. Thus Gortler vortices trapped in 
the wall layer have large wavelength with a = 0(l/(M 2 Mj)); whilst Gortler vortices trapped 
in the temperature adjustment layer have shorter wavelength and have 


a = 0(j r ). (3.8) 

It can be shown that the growth rates for these two modes are respectively of order one and 
order \J M 2 M \ . Thus the mode trapped in the temperature adjustment layer is the most 
dangerous mode and will be studied in the rest of this paper. 


4 The inviscid mode trapped in the temperature ad- 
justment layer 

As in the case when there is no pressure gradient (see Hall and Fu (1989)), an important 
term in the perturbation equations is that involving L(y) in the y-momentum equation. It 
is easy to show with the aid of the basic state properties listed in the Appendix that 

L(v) = + (4.1) 
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where the expression for B(X) is given in the Appendix. This term requires the Gortler 
number to be of order M 2 (the wall curvature should first of all overcome the effective 
negative curvature of the basic state if Gortler vortices were to exist at all). Hence we scale 
G by writing 

G=M 2 G. (4.2) 


With the order of G fixed, we can now use the basic flow properties to deduce the relative 
orders of the perturbation quantities. It is found that 

± = OisfMM,), V = 0(W) = 0(^0), P = 0(^-0), (4.3) 

and that U is of such an order that it does not enter the leading order analysis. The large 
growth rate implied by (4.3a) is driven by the term L(v) + | kGu 2 in the y-momentum 
equation. This makes the viscous terms in the perturbation equations unimportant and so 
the resulting eigenvalue problem is essentially inviscid. 

With the use of (4.1) and (4.2), we have 

L{v) + 1 -kGu 2 = Qk(*)G - B(x)j M 2 . 


Later we will show that if G is chosen such that the coefficient of M 2 in the above equation 
vanishes and if the crossflow is zero, the 0(\/M 2 Mi) growth rate will become zero and 
| k(x)G = B(x ) then defines the leading order neutral position (this is not true, however, 
when the crossflow is non-zero, as we shall show later). In the neighbourhood of this leading 
order neutral position, Gortler vortices have a smaller growth rate. Such near-neutral modes 
are dominated by viscous effects (see Hall and Fu (1989)). 

We now turn to the determination of the order of crossflow amplitude. It can be seen 
from (3.2)-(3.7) that this is achieved by investigating the terms involving the operator 

d d . _ 

u— + v— + taw. 
ox oy 

With the use of the asymptotic expressions listed in the Appendix, this operator can be 
written in the temperature adjustment layer as 


U,B l {\ + ■••)• (4.4) 

to leading order. It may appear that crossflow effects will be felt when djdX ~ aW\, but it 
will be shown later that the leading order constant term in w does not affect the leading order 
growth rate, see Hall (1985) for a similar result for the incompressible problem. Therefore, 
the appropriate size of the crossflow is fixed by letting 


d 

dX 


m 2P 


aWi 
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and in view of (3.8) and (4.3a), this gives 


W\ 




M 3 



+ 2 / 3 ' 


We thus introduce an 0(1) parameter 6 by 


Wi = 


N 




M 3 


Mi 




(4,5) 


to characterize the scaled crossflow amplitude. With the order of the crossflow fixed, we now 
look for the following form of WKB solutions for (3.2)-(3.7): 


V = 



(%(*,$) + ■•■)£ + ••■, 


w = \[j^( w o( x ,y) + ---)E + ..-, 

6 = (0 o {X,y) + ---)E + ---, (4.6) 

P = (Po(X,y) 4 ) E -f , 

where E , defined by 

£ = exp A (X)iX - iuW, f . (4.7) 

takes care of the fast variation of the perturbation quantities in the streamwise direction. 
Here A(X) is in general complex, the real part of it characterizes the growth rate of the 
disturbance. In each of the expansions in (4.6), the first ellipse represents the higher order 
terms which tends to zero as M — * oo, whilst the second ellipse represents the meanflow 
correction and higher harmonics. We note that the larger constant term iaW\ in (4.4) has 
been absorbed by inserting the second term in the expression for E. Thus, as remarked 
below (4.4), the term iaW\ in w only contribute a sinusoidal coefficient to the disturbance 
without affecting the growth rate. 

On substituting (4.6) into the linear perturbation equations (3.2) — (3.7) and make use of 
the basic state properties listed in the Appendix, we obtain, by keeping only those leading 
order terms, 

T x dV 0 _ 

“ ~T~d$ + lkW ° ~ °’ ( 4 - 8 ) 

C \ BP 

(A - ik8eP)V 0 - -0 O - = 0, (4.9) 

T p\ oy 
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(4.10) 


-Pa + S<e»Vi + £(A - .He*)VTo = 0, 

Pi T 2 T 

- Se>l' 0 + (A - .Me’OO.j = 0. 


(4.11) 


It can be seen that the function G defined above represents a balance between the positive 
wall curvature and the negative curvature due to the basic state. 

By eliminating W 0 ,9 0 and P 0 from (4.8) — (4.11) in favour of V 0 and introducing a new 
independent variable £ by 

C = e®, ’ (4.13) 


we obtain 


where 


,d 2 V 0 , C(l-C)dV c 


+ "(OK = 0 , 


l+c d( 


itS C(1 ~ C) , ,2 f GC n , na 

^ A - ikS( 1+C + {(A -ik6C) 2 ^ + ^ 


(4.14) 


(4.15) 


It is convenient to eliminate the first order derivative term from (4.14) by introducing V 
through 


1 + C 


Equation (4.14) then reduces to 


^ + {i + §TcF + H(c) } v = 0 ’ 


(4.16) 


where the ‘hat’ over V has been dropped to simplify notation. Under the transformation 
£ = e®, — oo < y < oo becomes 0 < ( < oo and therefore equation (4.16) is to be solved 
subject to the asymptotic boundary conditions V — * 0 as £ — > 0, oo. These conditions 
should be made more precise for numerical calculation and it can be deduced from (4.16) 
that the correct decay behaviour is 


v 


[1 + ^ + fc 2 (£ - 2)] + 0(C 2 )} as C - 0, 


2 k+ 1 


V — ► ( k e as £ — > oo. 


(4.17) 


The eigenvalue problem (4.16) and (4.17) determines the relationship between A (the real 
part of which is the growth rate), the scaled crossflow amplitude 6 and the scaled wavenumber 
k. In addition, G, defined by (4.12), can also be taken as a parameter for fixed X. We note 
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that for each set of fixed values of 6 , k and G, if (A, V) is an eigenvalue set of (4.16) then so 
is (—A*, V*), where an asterisk denotes complex conjugate. Therefore to each unstable mode 
there corresponds a stable mode. Also, for each set of fixed values of k and G , if (^, A, V) 
is an eigenvalue set of (4.16) then so is ( — S , — A,V). Hence we can restrict our attention to 
positive crossflow parameters and will do so for the remainder of this paper. 

Before we solve this eigenvalue problem for the non-zero crossflow cases, let us first note 
that when there is no crossflow, <5 = 0 and then (4.14) can be written as 


d 

d( 


C dV 0 \ k 2 
(i + C) 2 dCJ C(1 + C) 2 


-(i + 0 2 + 


§<}v„ = 0. 


(4.18) 


Equation (4.18) subject to Vo vanishing at ( = 0, oo is a standard Sturm- Liouville 
problem which can be taken to be an eigenvalue problem for G / A 2 for each fixed k. Its 
solution is of the form 

G . .o G 


Y* = F(k), i.e. A 2 = 


F(kY 


(4.19) 


where F{k) is a positive function of the scaled wavenumber k and can easily be determined 
by solving (4.18) numerically using the fourth order Runge-Kutta method. This solution 
will be discussed further in the following section. 

It follows from (4.19) that 


A 2 > 0 if G > 0, A 2 < 0 if G < 0. 

If we denote the location where G = 0 by X ni then A is real on one side of X n and is purely 
imaginary on the other side. This implies that as X = X n is crossed, the behaviour of the 
disturbance changes from being oscillatory to being exponentially growing, thus we may 
define X = X n as the neutral position and in view of the definitions (4.2) and (4.12b), the 
corresponding neutral Gortler number G then expands as 

G — 2i?(Xn ) m 2 _j_ J 1 ig] ier order terms, (4.20) 

K{X n ) 

where the higher order terms can be obtained by following the procedure given in Hall and 
Fu (1989). However, our purpose in the present paper is not to determine the expression 
for the neutral Gortler number, but instead we wish to investigate how the existence of 
a crossflow would affect the growth rate of the inviscid mode described by the eigenvalue 
problem (4.16). To be more specific, one of our aims in this paper is to determine how 
the growth rate— wavenumber relation shown in (4.19) for zero crossflow is modified as the 
crossflow is gradually increased. 
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5 Numerical results for non- neutral modes 


The eigenvalue problem (4.16) and (4.17) has been solved by using the fourth-order Runge- 
Kutta method to determine the dependence of the growth rate (the real part of A) on the 
wavenumber k with 6 fixed or on the crossflow amplitude S with k fixed. The boundary 
C = 0 is chosen to be some small number and ( = oo is replaced by a large number which is 
continuously adjusted for varying k and 6 values (since as we will see later the eigenfunctions 
become more and more concentrated around the singular point as we approach neutral 
stability). At each end, V is chosen arbitrarily and dVjdC^ is then calculated from (4.17). 
We shoot from the two ends and match the solutions at an appropriate middle point where 
we impose the condition that (1 !V)dV/d( > has the same valye from either side. 

For increasingly small wavenumbers, the eigenfunctions become more and more spread 
into the whole boundary layer (as we would expect since as k — » 0 we are approaching the 
wall layer mode) and in particular the eigenfunctions vary rapidly near the wall. In this case, 
it is found convenient to introduce a scaled variable £ through ( = e^ for the calculations in 
the interval to the left of the matching point. 

In order to be able to continue the calculations even when the mode becomes neutrally 
stable (i.e. when A becomes purely imaginary), we also introduce an indentation in our 
integration path. In order to see whether the indented path should be above or below the 
real axis, we first note from (4.15) that the complex singular point of the governing equation 
is at Q given by 

C = ^(Ai-V), (5.1) 

where A r + i\{ = A. For growing Gortler vortices A r > 0 and so this singular point lies below 
the real axis. Thus if the neutrally stable mode is taken to be approached from the unstable 
mode, then the integration path should be indented from above. In our calculations, the 
indentation is chosen to be a triangle and checks were duly made to make sure that it does 
not make any difference whether the indentation is on the left or right hand side of the 
matching point and how large the indentation is. 

We now present our numerical results for the non-neutral modes. All of the following 
results correspond to G = 2 and are for the first (i.e. the fundamental) mode. Other choices 
of positive G values can be shown to yield qualitatively similar results. The cases G < 0 
will be discussed in the final section. We note that although G < 0 is a possibility, cases 
corresponding to G > 0 are more important since when (5 = 0, positive values of G give rise 
to unstable Gortler vortices over the whole wavenumber regime, whilst if G < 0, no unstable 
Gortler vortices can exist. 

We shall first investigate the dependence of the growth rate on the wavenumber with 
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5 fixed. In Fig. 2(a) we have shown such growth rate curves for a number of 5 values, 
including the zero-cross flow case <5 = 0. We see that for 5 = 0, the growth rate tends to 
zero as k -> 0 and tends to a constant as k -+ oo. The former behaviour is required as the 
present mode should match with the wall layer mode as k — > 0 and the wall layer mode has 
smaller growth rate. We now show that the latter behaviour is governed by a WKB solution 
of the governing equation (4.18). 

When 5 = 0, Gortler vortices are governed by equation (4.18). As k — » oo, we look for 
WKB solutions of the form 

V 0 ~ exp (k j E{Qd(). 

Equation (4.18) requires that 

< ! E 2 (C) = )rC- (1 +<)’■ 

It then follows that in general there will be two first order turning points where E(() = 0. 
But the first mode which we are looking for must correspond to these two first order turning 
points coalescing to form a single second order turning point. This occurs when 

(f -2) ! -4 = 0, 

which yields the only acceptable result 


The corresponding second order turning point is then at £ = 1. In order to determine 
the higher order correction terms to the growth rate, we have to carry out the standard 
asymptotic analysis for second order turning point problems. Such an analysis shows that 
in the zero crossflow case A expands as 

A = rL - . k-' + 0(k-W). (5.2) 

z lb 

The dotted line shown in Fig. 2(a) corresponds to this asymptotic result. We note that the 
asymptotic analysis leading to (5.2) is very similar to that detailed in Fu, Hall and Blackaby 
(1992) for Sutherland’s law fluids and that the growth rate curves for the two viscosity laws 
have the same structure but they are fundamentally different from the growth rate curve for 
incompressible Gortler vortices. We also note that the growth rate curve for 5 = 0 given in 
Fig. 2(a) is also consistent with those growth rate curves, obtained by Dando and Seddougui 
(1992), which emerged as the Mach number is gradually increased. 
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As 8 is increased from zero gradually, we see from Fig. 2(a) that the growth rate is 
decreased at each fixed wavenumber and it no longer tends to a constant as k -* oo (this 
behaviour is also seen in Dando’s (1992) 0(1) Mach number calculations). Instead, we see for 
£ = 0.1 that there exists a critical wavenumber, k cr say, at which the growth rate vanishes. 
Such critical wavenumbers also exist for other 8 values (as can be seen from the exact solution 
to be discussed in the following section) although we are unable to show their existence in 
Fig. 2(a) due to numerical difficulty. For k > k CT , no solution can be found (since if a stable 
mode could be found, then an unstable mode could also exist). Thus at each 8 value, Gortler 
vortices are unstable in the wavenumber regime (0, k^), stable in the wavenumber regime 
(k CT ,o o) and are neutrally stable at k = k cr . 

Fig. 2(b) shows seven more neutral curves for the larger- 8 values indicated. We see that 
whilst for small 8 values crossflow effects are stabilizing over the whole wavenumber regime, 
this is no longer the case for large enough 5 values. Fig. 2(b) shows that crossflow effects are 
actually destabilizing over the small wavenumber regime whilst remaining stabilizing over 
the large wavenumber regime. We note that as 8 — ► oo, the critical wavenumber k cr at which 
A r = 0 tends to a definite limit. We shall show in the next section that this limiting value is 
given by 1/4. Hence we can deduce that for large 8 values, unstable Gortler vortices can only 
exist in the wavenumber regime (0, |) and that for any Gortler vortex with a wavenumber 
lying in the region (V,oo), there always exists a 8 value at which the vortex is completely 
stabilized. 

For Gortler vortices with wavenumber lying in the region (0, V), numerical calculations 
show that A = 0(8), as 8 — > oo. Therefore in the large crossflow limit, we can look for the 
following asymptotic solutions for (4.16) and (4.17): 


A — A 0 6 + Aj -f- A 2 8 1 + • ■ • , 


V(X, 0 = Vo(Jf, 0 + 1 V,(X, C) + • ■ • . (5.3) 

By substituting (5.3) into (4.16) and (4.17), it is easy to show that A 0 and Vo are determined 
by the eigenvalue problem 


.t &Vo 
K 3( 2 




ik 


4 (1 + () 2 A 0 - ik( 1 +C 


C(1 ^ -k 2 (\ +c) 2 W 0 = o, 


^^‘{’-TTSd i + f- 2 * 2 l + o ( c’)}, asC-0, 

Vo -> CV*<, as C -> oo. 


(5.4) 


We see that this reduced eigenvalue problem is now independent of G and hence the Gortler 
vortices in the wavenumber regime (0, |) can appropriately be identified as crossflow vortices 
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since in the large crossflow limit, their existence is dictated by the crossflow instead of by the 
wall curvature effects. The eigenvalue problem (5.4) can be solved numerically to determine 
the dependence of Ao on the wavenumber k, but we shall not present the corresponding 
results here since these results merely confirm the large crossflow behaviour of those growth 
rate curves shown in Fig. 2(b). 

In Figures 3(a, b), we have shown the dependence of the imaginary part of A on the 
wavenumber for the same 6 values as those used in Figures 2(a, b). We can see that A, > 0 
for A r > 0, which implies that when A r = 0, the critical point defined by (5.1) lies in the 
interval of integration and thus the neutral mode (to be studied in the following section) has 
a critical layer structure. 

In order to see the dependence of the growth rate on the crossflow for fixed wavenumbers, 
we have shown in Figures 4(a,b) the growth rate curves in the (A r , 8 ) plane for a selection of 
k values. We can immediately draw the following conclusions. For each fixed k value greater 
than 1/4, there exists an unique critical 8 value, 8 cr say, such that A r > 0 for 0 < 8 < <5^, 
A r = 0 at 8 = 8 cr , and there is no solution for 8 > 8 cr . For k < 1/4, we always have A r > 0, 
which confirms the findings shown in Fig. 2(b) that Gortler vortices with wavenumber lying 
in (0, j) are always unstable. For k = 1/4, neutral stability is attained at 8 cr = cc. Figures 
5(a,b) show the dependence of the imaginary part of A on the crossflow for the same k 
values used in Figures 4(a,b). Here again we have A; > 0. In Figures 6(a,b), we have shown 
the evolution of the real and imaginary parts of V for fixed values 8 — 0.2, G = 2 as k is 
increased. We can see that a critical layer structure is gradually emerging as k is increased 
towards the neutral value k cr . 


6 The neutrally stable mode 


Numerical results presented in the previous section suggest that to each 8 value, there cor- 
responds a critical k value at which neutral stability occurs and for each k value greater 
than 1/4, there again exists a critical 8 value at which neutral stability is achieved. In this 
section, we shall determine the relation between k and 8 when neutral stability occurs. 

When Gortler vortices are neutrally stable, we have A r = 0 and so we may write 


A — i A, 


( 6 . 1 ) 


where A is real and from our previous numerical results we know it is also positive. In the 
present neutral case, the governing equation (4.16) may be written as 


2 d 2 V 

d ( 2 


+ 



C(i-C) 
(1 + C) 2 


1 C(1-C) 

C - Cc 1+C 


cc 

S 2 (( - Cc) 2 


-* 2 (i + o 2 K = o, 


) 


( 6 . 2 ) 
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where 




(6.3) 


This equation has a regular singular point at ( = Cc and in the neighbourhood of this singular 
point the two Frobenius solutions can easily be shown to be given by 


V = V ± = ((-( c )^ ±l/) 1 + 


1 


-( 


1 -Cc 


G 


(1 ± U)( c 1 + Cc (c* 2 


)(C - Cc) + 


where 


v — 


\ 1 + Cc* 3 \ 


AG 


1 + 


AkG 

7T' 


(6.4) 


(6.5) 


In general, the solution in the neighbourhood of the singular point would be given by AV + + 
BV-i a linear combination of V+ and V-, but for the present neutral problem, we can show 
that either A or B must be zero. 

To the above end, we use Mile’s (1961) approach and consider a function r(X,() defined 

b y 

r(X,C) = (V'n,, (6-6) 


where a subscript i 

signifies taking the imaginary part and a prime denotes partial differen- 

_ 

tiation with respect to C- It is straightforward to show with the aid of ( 6 . 2 ) that dr/d( — 0. 
Thus r can only be a function of X only. But we know t(X, 0) = t ( oo ) = 0 and so 


* 

</■% 

III 

o 




However, if 



- 


V = AV+ + BV-, 

(6.7) 

- 

we then have 



7 


T = (AB*Vl VI + A*BViV*)i. 

(6.8) 

— 

As we have argued 

in the previous section, if the neutral mode is 

taken to be the limit of 

- 

the unstable modes, then the branch cut should be in the + plane. 

We therefore have 

= 


* if / 0 C - Cc > 0, 

Ar g(C Cc)-j x? C-Cc<0, 


_ 

and 





(C - Q* xiv) = 1C - Cc|^ (1±l/) e*^ (1±i/)5 , 

(6.9) 

- 

where 

/ o, C - Cc > o, 

1 l, C - Cc < o. 

(6.10) 

- 

Hence we have 



— 


V* = V±e“ i,r(1:l:l,)S . 

(6.11) 

- 
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On substituting (6.1 1) into (6.8), and making use of the result V+V- — V + V!_ = //, we obtain 

r = {AB*i/e ,,r(1+t/)5 }.. 

Hence by evaluating the above expression on the left and right hand sides of £ c , we obtain 

{AB*)i = 0 and (AB*e iir ^ + % = 0, 

which imply that AB* = 0 whenever sin7r/^ ^ 0 and so either A or B must vanish if v is not 
an integer (our numerical results show this is also true even when v is not an integer). 

Once we have established this result, we can determine whether the solution is AV+ or 
BV- with the aid of the numerical solutions for the neutral mode (as the limit of our non- 
neutral calculations discussed in the previous section). The numerical results for G > 0 show 
that the eigenfunctions are singular at the critical point and thus we deduce from (6.4) that 
the solution is always given by BV 1 (since v > 1 when G > 0 and V+ is regular). This can 
also be verified by calculating the phase jumps across the critical point of the eigenfunctions; 
we then find that the phase jump is always tt(1 — ^)/2 instead of 7r(l + //)/2. Here v is 
calculated from (6.5). In the course of calculating these phase jumps for a selection of k 
values, we noticed the possibility v = 4k and this was then verified by calculations for a wide 
range of wavenumbers. This result then points to the possible existence of an exact solution 
to the neutral eigenvalue problem (6.2). After some manipulation, such an analytical solution 
is found to indeed exist and is given by 

is = 4k, Cc = -p 6 = 4A 

16 Gk = (16 k 2 - 1 )8\ (6.12) 

i/ = fichc-Cc)^- u (c+i)" 1 e- fcC . 

We note that each component in the product on the right hand side of (6.12c) has a distinct 
role: they in turn represent the behaviour of the eigenfunction at the four regular singular 
points ( = 0, ( c , —1, and oo. 

The exact solution (6.12) can explain quite a few points about our numerical results. 
Firstly, we can see from (6.12d) that k = 1/4 is an important value. For any fixed G > 0, 
neutral states can only be obtained for k > 1/4 since (6.12d) does not have a solution for 
8 otherwise. When k = 1/4, the neutral state is obtained at 8 — oo. Secondly, for a given 
8 value, there exists a neutral k value (i.e. a k cr at which A r = 0) and for a given k value 
greater than 1/4, there exists a neutral 8 value. All these conclusions certainly agree with 
the numerical results presented in Figures 2-5. Finally, we see from (6.12b) that even for 
moderate values of k, ( c can be quite small and in these cases the eigenfunctions become 
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trapped near the wall, which give rise to some numerical difficulties. This is why we have 
not extended the growth rate curves to the neutral points in Fig. 2(a) although we have 
done so in Fig. 2(b). 

After having succeeded in obtaining the above exact neutral mode solution, we also 
considered the neutral mode eigenvalue problem formulated, but unsolved, by Bassom and 
Hall (1991). As remarked by these authors, their equation (2.14b) is just a scaled version of 
their equation (4.10c) and it governs the large wavenumber inviscid neutral mode. It suffices 
to consider their equation (4.10c), which can be written as 


where tp c = —0/ A. Using the same procedure as that used in obtaining (6.12), we found 
that (6.13) also has an exact solution, given by 

A = -/3 = -L, t? 0 = AW-jJ-'e-*, (6.14) 

where A is an arbitrary constant. 

The eigenvalue problem associated with (6.13) arises in Bassom and Hall (1991) where 
the re-emergence of the Gortler vortex mode is considered at high scaled crossflow and 
vortex wavenumbers. In that paper it was found that at sufficiently high crossflow values 
the Gortler mode is stabilized over a finite band of nonzero wavenumbers. This stable 
band of wavenumbers increases with the crossflow; its left-hand limit approaches that of the 
stationary crossflow vortex of Gregory, Stuart and Walker (1955), while its right-hand limit 
is described by the inviscid eigenvalue problem associated with (2.14b) of Bassom and Hall 
(1991). The result (6.13) shows that, in the notation of Bassom and Hall (1991), the upper 
neutral wavenumber is related to the crossflow parameter A by 

A = 2V2 a I/2 ( /r I/ 2 G ,/2 . (6.15) 

V w (°)/ 

Here R and G are the Reynolds and Gortler numbers whereas the unperturbed boundary 
layer is (u, R~ l t 2 v, \w). The prediction (6.15) is in excellent agreement with the numerical 
results of Bassom and Hall (1991) thus for example the result (6.15) differs from the computed 
value of a by less than 5% when A R l ^ 2 G~ 1 ^ 2 = 4.75. 


_ (1 + - + „ ! 

dtp 2 \ xp(xp - xpc) A 2 1p(lp - xpc) 2 
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7 Conclusion 


So far, we have consistently assumed that G > 0. We have shown that if 8 = 0, unstable 
Gortler vortices can exist in the whole wavenumber regime. But as 5 is increased gradually 
from zero, this (0, oo) unstable wavenumber regime contracts to (0, k ) where k is related to 
8 by (6. 1 2d) and tends to 1/4 as 8 tends to infinity. Any Gortler vortex with a wavenumber 
greater than 1/4 could be completely stabilized by a large enough crossflow in the basic 
state. We have also shown that Gortler vortices over the wavenumber regime (0, j) are 
always unstable and for large crossflows they are actually destabilized by the crossflow. 
Since such Gortler vortices for large crossflows are no longer controlled by wall curvature 
effects, we have identified them more appropriately as crossflow vortices. 

To complete our discussions, we now investigate what happens if G < 0 . In the zero- 
crossflow case, result (4.19) shows that if G < 0 , the inviscid mode has \ r = 0 and so the 
growth rate curve in the (k,\ r ) plane is simply the £-axis. We now investigate how this 
curve (i.e. the k-axis) evolves as we increases 8 gradually. We note that the neutral solution 
(6.12) is also valid for G < 0. When G < 0, we deduce from (6.1 2d) that for each 8 value, 
there exists an unique k , fc CT say, at which A r = 0, and we must have k cr < 1/4. This is 
verified by the numerical results shown in Figure 7(a,b) where we have shown two typical 
growth rate curves corresponding to 8 = 5, 7 and to G = —2. It can be seen that for 8 > 0 
and G < 0, the only unstable wavenumber regime is now (0, k cr ) with k CT < 1/4, which is 
even narrower than that for G > 0. This is a result of the crossflow trying to destabilize 
and the negative resultcmt curvature trying to stabilize these (crossflow) vortices. It can be 
deduced from (6.1 2d) , and it is also confirmed by Fig. 7, that k cr tends to 1/4 from below as 
8 — ► oo, in contrast to the situation for G > 0 where k cr tends to 1/4 from above as 8 — » oo. 

Figures 8(a,b) show that a crossflow vortex at a fixed 8 value and with wavenumber 
lying in the region (0, |) can be completely stabilized by decreasing G , and so although the 
crossflow vortices in the wavenumber regime (0, A) are always unstable for G > 0, they can be 
nevertheless completely destroyed by a large enough resultant negative curvature (represented 
by G). Therefore, crossflows stabilize the unstable Gortler vortices with wavenumber lying 
in the region (1/4, oo) and the negative resultant curvature stabilizes the unstable crossflow 
vortices with wavenumber lying the region (0, 1 /4). 

To summarize, we have shown in this paper that in the hypersonic context, the crossflow 
which can completely stabilize the inviscid Gortler vortices considered must necessarily be 


of order 


/T 1/2 


M 3 

M2+ 2 ^ 


in terms of the original physical variable. Within this order, we have shown at precisely what 
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crossflow value a Gortler vortex with a given wavenumber will be completely stabilized. It 
is seen that large wavenumber Gortler vortices can be stabilized more easily than smaller 
wavenumber Gortler vortices. This may imply that in practical situations small wavenumber 
Gortler vortices could be excited more easily than Gortler vortices with larger wavenumbers. 

It is also seen that crossflow effects on Gortler vortices are fundamentally different for 
incompressible and hypersonic flows. For incompressible flows, Bassom and Hall (1991) 
showed that only at sufficiently large crossflow values, were Gortler vortices stabilized over 
a finite band of wavenumbers, whilst we have shown in this paper that for any non-zero 
crossflow, hypersonic Gortler vortices are stabilized over a semi-infinite wavenumber regime 
(k CT} oo), where k cr is determined by the exact solution (6.12d) (with k cr identified with the 
k there). 

We have also shown that in order to quantify vortices’ response to crossflow effects, the 
scaled wavenumber regime may be neatly divided into two regions: (0, 1/4) and (1/4, oo). 
We have identified the vortices wavenumber lying in the first region as crossflow vortices 
since in the large crossflow limit, their existence is no longer supported by curvature effects 
but instead by the crossflow. We have also found the exact solution to the neutrally stable 
mode. This exact solution should be useful for studying the nonlinear development of Gortler 
vortices which is now under consideration. 
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Appendix: Basic state properties in the temperature adjustment 

layer 
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Figure 1 . A typical example of a three-dimensional boundary layer: a boundary layer over 
a swept aerofoil. 
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Fig. 5(a) 



Fig. 5(b) 



Figure 5. Dependence of the imaginary part of A on the crossflow amplitude for G — 2. (a) 
For k = 0.5, 0.6, 0.7,..., 1.2; (b )for k=0.1, 0.15, 0.2, 0.25, 0.3. 
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Fig. 6(b) Eigenfunctions for 6=0.2, G=1 



Fig. 6(a) Eigenfunctions for 6 = 0.2, G=1 



0.0 0.2 0.4 0.6 0.8 1.0 


Figure 6. Variation of the eigenfunction for G = 2 and 8 = 0.2 as k is increased towards 
the neutral value, showing the emergence of a critical layer structure, (a) The real part of 
V; (b) the imaginary part of V. 
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